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PREDICTION  OF  CLASSICAL  FRACTURE  INITIATION  TOUGHNESS 
INTRODUCTION 

Linear  elastic  fracture  theories  have  been  used  successfully  to  predict  brittle  fractures.  How¬ 
ever,  brittle  failures  represent  only  a  small  portion  of  failures  due  to  the  presence  of  a  crack.  High 
toughness  materials  may  exhibit  large  amounts  of  local  deformation  prior  to  fracture.  In  order  to 
accurately  predict  failure  in  high  toughness  materials,  damage  and  fracture  theories  which  address 
ductile  behavior  are  required. 

Many  failure  criteria  have  been  proposed  for  ductile  materials.  Predominant  among  the  failure 
criteria  proposed  are  contour  integrals  such  as  J  [1]  and  crack  tip  opening  displacement  [2].  Some 
early  work  in  the  determination  of  ductile  failure  criteria  consisted  of  use  of  correction  factors  to 
extend  the  region  of  applicability  of  linear  elastic  fracture  mechanics  by  representing  the  effects  of 
plasticity  around  the  crack  tip  [3,4],  However,  only  a  limited  amount  of  plasticity  can  be  approxi¬ 
mated  by  such  correction  factors.  In  other  instances  linear  elastic  material  damage  and  fracture  cri¬ 
terion  evolved  into  ductile  criterion  through  enhancements  which  broadened  the  range  of  applicability 
of  the  failure  criterion.  The  resistance  curve  analysis  method  is  one  such  evolved  damage  and  frac¬ 
ture  criterion.  Initial  resistance  curve  analyses  were  based  on  a  plasticity  corrected  stress  intensity  fac¬ 
tor  [5].  Resistance  curve  concepts  were  then  extended  to  the  case  of  limited  plasticity  and  perfecdy 
plastic  material  response  with  the  incoporation  of  J  as  the  failure  criterion  as  shown  permissible  by 
Landes  and  Begley  [6].  Additional  work  concentrated  on  developing  J  for  conditions  of  full  plasticity 
extending  for  a  relatively  large  area  around  the  crack  tip  [7,8].  With  each  stage  of  its  development, 
the  resistance  curve  concept  has  been  extended  in  an  effort  to  account  for  increased  material  ductility. 

In  general,  when  there  is  a  relatively  small  amount  of  plasticity  around  the  crack  tip,  failure  cri¬ 
teria  based  on  linear  elastic  concepts  can  be  used  successfully.  However,  in  instances  where  large 
amounts  of  plasticity  occur  prior  to  failure,  analysis  methods  based  on  linear  elastic  concepts  are  no 
longer  capable  of  approximating  the  state  of  stress  and  strain  near  the  crack  tip.  In  order  to  accurately 
represent  the  state  of  stress  and  strain  near  the  crack  tip  the  analysis  method  and  failure  criteria  used 
should  explicitly  include  large  deformation  constitutive  response. 

The  explicit  consideration  of  ductile  material  behavior  prior  to  fracture  does  not  in  and  of  itself 
alleviate  all  limitations  which  may  be  associated  with  a  particular  fracture  mechanics  analysis  tech¬ 
nique.  Present  day  fracture  mechanics  analysis  techniques  have  limitations  which  must  be  taken  into 
consideration  [9].  Limitations  common  to  fracture  mechanics  analysis  techniques  include  the 
mathematical  prediction  of  physically  unattainable  conditions  such  as  the  existence  of  a  stress,  strain 
and  strain  energy  density  singularity,  correct  material  characterization  and  the  appropriateness  of  the 
chosen  fracture  criterion.  The  material  characterization  used  for  ductile  fracture  should  include  an 
accurate  description  of  large  strain  response.  The  failure  criterion  for  ductile  fracture  should  not  be 
limited  by  the  amount  of  local  deformation  which  occurs  prior  to  fracture.  Other  areas  of  concern 
which  may  exist  for  specific  problems  include  availablity  of  suitable  data  to  adequately  characterize 
material  deformation,  effects  of  loading  rate  and  history,  interpretations  of  material  data  scatter  and 
defect  characterization. 

In  addition  to  the  various  limitations  listed  above,  a  feature  common  to  many  fracture  criteria  is 
the  assumption  that  failure  occurs  at  a  dominant  crack  or  defect  which  is  observed  or  assumed  to  exist 
in  the  structure.  Since  failure  is  assumed  to  occur  at  a  crack  or  a  defect,  these  features  must  be 
present  in  the  test  specimen  used  to  determine  the  critical  values  of  the  parameters  in  the  fracture  cri¬ 
terion. 
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A  material  damage  and  failure  criterion  which  explicitly  addresses  ductile  behavior  can  be 
developed  based  on  a  general  continuum  approach  which  utilizes  the  constitutive  relations  to  describe 
deformation,  damage  and  fracture  by  the  continuous  point  variables  of  stress,  strain  and  energy  den¬ 
sity.  Continuum  material  toughness  [10]  is  one  such  material  damage  and  fracture  criterion  which  is 
based  on  a  general  continuum  approach  and  is  the  criterion  used  in  the  current  analysis.  For  a  struc¬ 
tural  analysis,  the  geometry  and  boundary  conditions  which  define  the  problem  and  the  constitutive 
relations  representing  the  material  are  used  in  a  computational  solution  for  the  stress,  strain  and 
energy  density  fields.  The  critical  locations  in  the  structure  are  identified  from  the  computational 
simulation  and  the  load  or  time  to  failure  is  determined  from  the  local  stress,  strain  and  energy  den¬ 
sity  histories.  No  a  priori  assumption  as  to  the  location  or  mode  of  failure  is  required.  The  failure  cri¬ 
terion  used  is  the  strain  energy  density  required  to  produce  material  fracture  at  the  continuum  scale. 
The  use  of  an  energy  density  as  a  failure  criterion  is  consistent  with  the  concepts  presented  by  Freu- 
denthal  [11]  on  material  behavior  and  scaling  considerations.  Previous  work  by  Gillemot  [12]  used 
analytical  and  empirical  techniques  to  determine  the  strain  energy  density  per  unit  volume  at  fracture 
for  cylindrical  tensile  specimens.  The  strain  energy  per  unit  volume  absorbed  by  the  material  up  to 
the  instant  of  fracture  was  calculated  from  global  specimen  response  and  deformation  geometry.  In 
separate  studies  Matic  et  al  [13]  have  used  computational  simulation  along  with  the  experimentally 
determined  point  of  failure  to  determine  the  critical  strain  energy  density.  Calculations  of  critical 
strain  energy  density  values  from  computational  simulations  are  able  to  address  large  deformation 
continuum  toughness  in  a  manner  appropriate  for  structural  integrity  analyses.  Applications  to  struc¬ 
tural  integrity  prediction  which  demonstrate  the  feasiblity  of  continuum  material  toughness  as  a  failure 
criterion  include  investigations  by  Matic  and  Joiles  [14,15]  into  the  defect  tolerance  of  welded  com¬ 
ponents  and  the  role  of  material  geometry  and  applied  load  on  weld  system  performance. 

In  order  to  demonstrate  the  fundamental  nature  of  continuum  material  fracture  toughness  as  a 
material  damage  and  fracture  criterion,  a  computational  simulation  of  ductile  fracture  initiation  in  a 
compact  specimen  of  HY-100  steel  under  monotonic  loading  is  performed.  Results  from  a  fracture 
initiation  experiment  are  available  for  comparison.  Continuum  material  toughness  concepts  are 
reviewed  in  order  to  clarify  for  the  reader  the  choice  and  implementation  of  material  fracture  tough¬ 
ness  at  the  continuum  scale  as  a  failure  criterion.  The  critical  strain  energy  density  required  for  frac¬ 
ture  at  the  continuum  scale  and  the  scale  invariant  continuum  stress-strain  relation  are  determined  by 
appropriate  analysis  of  data  obtained  from  a  series  of  cylindrical  tensile  specimens  of  the  same 
material  as  the  compact  specimen.  The  continuum  stress-strain  relation  is  determined  from  hybrid 
computational-experimental  analyses  which  include  large  strain  response.  Nonlinear  finite  element 
methods  are  used  for  the  computational  analyses. 

CONTINUUM  MATERIAL  TOUGHNESS  CONCEPTS 


A  continuum  volume  of  material  undergoing  deformation  is  characterized  by  its  multiaxial  stress 
and  strain  state,  and  ,  respectively.  For  ductile  metals,  and  other  materials  which  exhibit  inelas¬ 
tic  deformation,  the  strain  state  is  a  function  of  the  stress  history  as  well  as  the  stress  state.  The  state 
of  stress  and  strain  may  be  related  by  a  constitutive  formulation  which  defines  the  strain  increment 
dfij  from  the  current  stress  state. 

The  strain  energy  per  unit  mass  at  a  given  instant  during  deformation  is: 
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where  p  is  the  mass  density.  The  energy  density  incorporates  both  stress  and  strain  into  a  fundamental 
quantity  relevant  to  thermodynamic  description  of  material  deformation  and  damage.  The  energy  den¬ 
sity  is  a  scalar  quantity  which  takes  into  account  all  components  of  the  stress  and  strain  tensors  in  a 
physically  consistent  manner.  Failure  of  the  material,  at  the  continuum  scale,  can  be  associated  with 
the  value  of  the  energy  density  at  which  fracture  occurs.  Thus,  the  material  toughness  may  be  defined 
as: 
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where  wr  is  the  critical  strain  energy  density.  The  value  of  wr  can  be  considered  as  a  material  pro¬ 
perty  for  engineering  analyses. 


For  ductile  metals,  the  material  density  varies  only  slightly,  even  over  large  deformations.  For 
this  reason,  it  is  possible  to  define  an  energy  per  unit  volume  density  as: 
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with  an  associated  critical  value 
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The  energy  per  unit  mass  is  fundamental,  however  the  energy  per  unit  volume  is  equally  appropriate 
for  constant  volume  deformation  processes. 


For  the  case  of  an  uniaxial  stress-strain  curve,  corresponding  to  a  one-dimensional  state  of 
deformation,  the  critical  strain  energy  density  corresponds  to  the  area  under  the  uniaxial  stress-strain 
curve: 


wc  =  $0  ade  (?) 

This  representation  is  desirable  for  use  with  traditional  constitutive  formulations  which  rely  on 
equivalent  uniaxial  stress-strain  curves.  In  the  general  case  of  multiaxial  stress  and  strain,  the  strain 
energy  density  is  the  sum  of  the  integrals  of  the  individual  stress  and  strain  components. 

MATERIAL  CHARACTERIZATION 

In  order  to  obtain  an  accurate  computational  simulation  of  fracture  initiation,  the  constitutive  for¬ 
mulation  used  must  be  valid  for  the  entire  deformation  range  between  zero  deformation  and  fracture 
initiation.  A  hybrid  computational-experimental  procedure  is  used  to  determine  the  constitutive  rela¬ 
tion  for  HY1-00  steel  [13].  Data  from  tensile  tests  performed  using  cylindrical  specimens  is  used  to 
characterize  the  material  behavior  of  HY-100  steel.  The  continuum  stress-strain  relation  determined 
from  the  tensile  tests  is  based  on  the  entire  range  of  deformation  observed  in  the  ductile  material. 
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This  hybrid  computational-experimental  approach  uncouples  the  influence  of  test  specimen 
geometry  from  the  observed  deformation.  This  results  in  a  unique  scale  invariant  continuum  stress- 
strain  curve.  The  solution  curve  for  HY-100  steel  is  shown  in  Figure  1  and  is  the  stress-strain  relation 
used  in  the  current  analysis. 

The  critical  strain  energy  density  as  defined  by  equation  (6)  or  (7)  is  determined  from  the  solu¬ 
tion  continuum  stressstrain  curve,  the  computationally  predicted  stress-strain  histories  at  the  point  of 
local  fracture  and  the  experimentally  observed  load-displacement  point  corresponding  to  tensile  speci¬ 
men  fracture.  Low  and  high  aspect  ratio  tensile  specimens  result  in  different  critical  strain  energy 
density  values  [13].  The  low  aspect  ratio  tensile  specimens  exhibited  multiaxial  deformation 
throughout  the  local  stress-strain  history.  In  contrast,  the  high  aspect  ratio  tensile  specimens  sustained 
a  predominantly  uniaxial  state  of  stress  prior  to  specimen  necking  followed  by  a  transition  to  multiax¬ 
ial  stress  and  strain  as  the  neck  develops.  The  strain  energy  density  at  the  point  of  material  fracture  in 
the  low  aspect  ratio  specimens  is  used  as  the  critical  strain  energy  density  for  the  current  analysis 
since  these  specimens  are  more  representative  of  the  multiaxial  stress  and  strain  history  ahead  of  the 
crack  tip  in  the  compact  specimen.  The  critical  strain  energy  density  coinciding  with  fracture  obtained 
from  the  low  aspect  ratio  tensile  specimens  for  HY-100  steel  is  882.6  MN-m/m3  (129,000  lb-in/in3). 

FRACTURE  INITIATION  EXPERIMENT 

A  compact  specimen  was  used  to  experimentally  determine  classical  measures  of  fracture  initia¬ 
tion  toughness  of  HY-100  steel.  Nominal  specimen  dimensions,  as  shown  in  Figure  2,  are  B  =  22.9 
mm  (0.90  in),  W  =  50.8  mm  (2.0  in)  and  a/W  =0.5;  where  B,  W  and  a  are  the  specimen  thickness, 
width  and  crack  length,  respectively.  A  nominal  2.54  mm  (0.10  in)  crack  measured  from  the  notch  tip 
was  produced  by  fatigue  prior  to  testing.  This  was  accomplished  with  a  servohydraulic  testing 
machine  in  load  control  mode  with  the  load  cycling  between  1.56  and  15.57  kN  (350  and  3500  lb). 
Prior  to  precracking  the  alignment  of  the  machine  was  checked  and  found  to  be  within  0.0508  mm 
(0.002  in).  After  the  precrack  started  to  grow  the  specimen  was  turned  over  as  necessary  to  keep  the 
front  and  back  face  crack  lengths  nearly  equal.  The  crack  lengths  at  the  front  and  back  face  were 
monitored  and  measured  with  traveling  microscopes. 

An  alternating  current  (AC)  electric  potential  measurement  technique  was  used  to  determine 
when  crack  initiation  occured.  The  placement  of  the  leads  on  the  specimen  are  shown  in  Figure  2.  A 
lock-in  amplifier  provided  a  2  Amp,  5  kHz  reference  signal  to  leads  A  and  measured  the  potential 
across  the  crack  mouth  at  leads  B. 

A  series  of  tests  were  performed  on  HY-80  steel  to  assess  the  reliablity  of  the  AC  electric  poten¬ 
tial  method  in  determining  fracture  initiation  in  ductile  steels.  These  tests  demonstrated  that  initiation 
of  crack  growth  is  coincident  with  attaining  a  minimum  in  the  measured  potential. 

The  specimen  was  loaded  under  displacement  control  using  a  servohydraulic  test  machine  at  a 
displacement  rate  of  7.1  x  10-3  mm/sec  (2.8  x  10-4  in/sec).  A  microcomputer  started  the  test,  per¬ 
formed  data  acquistion  of  load  and  potential  across  the  crack  mouth  and  also  controlled  a  video  digi¬ 
tizing  board  and  video  camera.  The  microcomputer  acquired  data  at  ten  second  intervals.  For  each 
data  set  a  magnified  image  of  the  crack  mouth  was  digitized  and  displayed  on  a  video  monitor.  The 
images  were  recorded  by  a  video  cassette  recorder  for  subsequent  measurements  of  the  crack  mouth 
opening  displacement  (CMOD).  As  soon  as  it  was  realized  that  the  potential  reached  a  minimum  the 
specimen  was  unloaded,  removed  from  the  test  machine  and  marked  by  heat  tinting  at  500 °F  for  90 
minutes.  The  specimen  was  then  cooled  in  liquid  nitrogen,  placed  in  the  test  machine  and  rapidly  bro¬ 
ken.  Examination  of  the  tinted  surface  indicated  a  negligible  amount  of  crack  growth.  Post  test  meas¬ 
urements  of  the  CMOD  were  performed  by  manually  measuring  the  distance  between  the  top  and  bot¬ 
tom  comers  of  the  crack  mouth  of  the  digitized  image  with  a  cursor.  The  resolution  of  this  technique 
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was  0.0404  mm  (0.0016  in).  The  instrumentation  is  shown  schematically  in  Figure  3.  The  experi¬ 
mentally  determined  fracture  initation  load  is  78.06  kN  (17,550  lb)  and  the  measured  CMOD  at  frac¬ 
ture  initiation  is  1.17  mm  (0.0459  in). 

FRACTURE  INITIATION  SIMULATION 

Fracture  initiation  in  a  compact  specimen  of  HY-100  steel  is  computationally  simulated  using  the 
ABAQUS  finite  element  code  [16]  and  the  concepts  described  earlier.  Nominal  test  specimen  dimen¬ 
sions  are  used  in  the  computational  simulation.  The  crack  in  the  compact  specimen  is  assumed  to  have 
a  straight  front. 

Two  and  three-dimensional  nonlinear  analyses  were  performed.  The  three-dimensional  analysis 
was  performed  when  results  of  the  two-dimensional  analysis  indicated  that  a  two-dimensional  solution 
was  inadequate  to  accurately  reflect  the  physically  observed  deformation. 

Two-Dimensional  Simulation 

The  two-dimensional  model  is  shown  in  Figure  4.  A  total  of  2534  C1D2,  CPE8  and  CPE8H  ele¬ 
ments  are  used  to  model  half  of  the  compact  specimen.  C1D2  elements  are  2-node  beam  elements 
used  for  applying  load  to  the  model.  CPE8  and  CPE8H  elements  are  both  8-node  continuum  ele¬ 
ments  with  quadratic  displacement  interpolation.  In  addition  to  quadratic  displacement  interpola¬ 
tion  the  CPE8H  elements  incorporate  an  independently  interpolated  linear  hydrostatic  stress  [17]. 
Lagrange  multipliers  are  used  to  couple  the  hydrostatic  stress  to  the  constitutive  relationship.  The  use 
of  such  hybrid  elements  in  regions  of  large  deformation,  such  as  near  the  crack  tip,  prevent  physically 
unrealistic  nodal  displacement  constraints  from  propagating  through  the  mesh.  In  a  model  comprised 
of  only  nonhybird  elements  these  unrealistic  displacements  artificially  increase  the  stiffness  for 
incompressible  or  nearly  incompressible  material  deformation. 

The  finite  element  mesh  is  more  refined  in  the  region  surrounding  the  crack  tip  than  in  the 
remaining  portion  of  the  model.  Each  element  along  the  crack  line  is  0.4  percent  of  the  total  crack 
length.  No  crack  tip  singularity  is  included  in  the  finite  element  model. 

Fifteen  C1D2  beam  elements  are  used  to  apply  load  to  the  compact  specimen  model.  The  beam 
elements  are  arranged  in  a  finite  element  model  is  loaded  by  application  of  a  point  force  at  the  center 
node  of  the  ’wagon- wheel’  loading  frame.  The  beam  elements  are  rigid  elastic  connectors.  Lateral 
translation  of  the  finite  elment  model  is  prevented  by  restraining  lateral  movement  of  the  center  node 
of  the  beam  element  loading  frame. 

A  nonlinear  static  plane  strain  analysis  was  performed  using  ABAQUS.  Both  geometric  and 
material  nonlinearities  are  included  in  the  analysis.  Geometric  nonlinearity  is  included  for  accurate 
evaluation  of  large  strains  and  rotations  which  may  occur  locally  prior  to  failure.  Use  of  Rik’s  algo¬ 
rithm  assures  numerical  convergence  of  the  finite  element  solution  independent  of  structural  stability. 
In  the  modified  Rik’s  algorithm  available  in  ABAQUS,  the  magnitude  of  the  applied  load  is  an 
unkown  to  be  determined  as  part  of  the  solution.  The  program  determines  increments  of  load  based 
on  increments  of  path  length  along  the  load-displacement  solution  curve  for  the  structure.  The  path 
length  increments  are  chosen  by  the  program  so  that  a  convergent  numerical  solution  is  obtained. 

The  finite  element  model  was  subjected  to  monotonically  increasing  load  until  fracture  initiation. 
Fracture  initiation  occurs  when  the  strain  energy  density  at  the  crack  tip  reaches  the  local  material 
toughness  determined  from  the  tensile  specimens  as  described  earlier. 
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The  strain  energy  density  profile  ahead  of  the  crack  tip  is  shown  in  Figure  5.  The  strain  energy 
density  is  observed  to  monotonically  decrease  as  the  distance  from  the  crack  tip  increases. 

The  global  load-displacement  curves  for  the  two-dimensional  finite  element  model  and  the 
experiment  are  shown  in  Figure  6.  The  displacement  shown  is  total  CMOD.  At  low  load  levels  there 
is  good  agreement  between  the  experimental  and  predicted  loaddisplacement  curves.  However,  at 
approximately  80  percent  of  the  fracture  initiation  load  the  experimental  Jesuits  indicate  a  reduction  in 
stiffness  which  is  not  observed  in  the  two-dimensional  finite  element  analysis  >esults.  The  difference 
in  stiffnesses  near  fracture  initiation  is  considered  to  be  a  result  of  two-dimensional  modeling.  Physi¬ 
cally  occuring  out  of  plane  displacements  and  their  effects  on  the  stress  and  strain  state  within  the 
compact  specimen  cannot  be  accurately  modeled  by  a  two-dimensional  plane  strain  analysis. 

Three-Dimensional  Simulation 

A  three-dimensional  finite  element  analysis  was  performed  in  order  to  obtain  a  more  accurate 
solution  of  the  global  response  of  the  compact  specimen. 

The  three-dimensional  finite  element  model  is  shown  in  Figure  7.  Symmetry  conditions  allow 
for  modeling  one-fourth  of  the  actual  specimen.  A  total  of  301  elements  were  used  to  generate  the 
model.  Even  though  fewer  elements  are  used  in  the  three-dimensional  finite  element  model,  the  two 
models  require  approximately  the  same  amount  of  computer  memory. 

The  three-dimensional  model  is  made  of  C3D20  and  C3D20H  elements.  These  are  20-node  con¬ 
tinuum  brick  elements  with  quadratic  displacement  interpolation.  In  addition  to  quadratic  displacement 
interpolation  the  C3D20H  elements  incorporate  an  independently  interpolated  linear  hydrostatic  stress. 
As  in  the  two-dimensional  analysis,  use  of  such  hybrid  elements  in  regions  of  large  deformation 
prevents  artificial  stiffening  of  the  finite  element  model. 

The  three-dimensional  finite  element  model  is  considerably  less  refined  than  the  two-dimensional 
model.  Each  element  along  the  crack  line  is  approximately  9.0  percent  of  the  crack  length.  The 
three-dimensional  model  has  local  refinement  at  the  crack  tip  and  near  the  outer  surface  of  the  com¬ 
pact  specimen.  There  are  seven  elements  across  the  thickness. 

A  nonlinear  static  analysis  was  performed  using  ABAQUS.  The  model  was  subjected  to  monton- 
ically  increasing  displacement  until  fracture  initiation.  As  in  the  two-dimensional  analysis,  geometric 
and  material  nonlinearities  are  included,  Rik’s  algorithm  is  used  to  guarantee  the  numerical  conver¬ 
gence  of  the  finite  element  solution  and  no  crack  singularity  is  included  in  the  finite  element  model. 

The  strain  energy  density  profile  across  the  specimen  thickness  at  the  plane  of  the  crack  is 
shown  in  Figure  8.  The  maximum  strain  energy  density  observed  in  the  finite  element  model  occurs 
at  the  midplane  of  the  specimen.  Therefore,  fracture  initiation  is  predicted  to  begin  at  the  midplane  of 
the  specimen.  Even  though  the  initial  crack  front  is  straight,  crack  initiation  as  represented  by  the 
strain  energy  density  is  nonuniform  across  the  specimen  thickness.  The  shape  of  the  strain  energy 
density  profile  across  the  specimen  thickness  is  reminiscent  of  the  curved  crack  front  observed  in 
compact  specimens  which  exhibit  ductile  tearing  as  shown  in  Figure  9.  The  three-dimensional  finite 
element  analysis  load-displacement  solution  is  shown  in  Figure  10. 

This  three-dimensional  finite  element  model  does  not  have  the  same  degree  of  refinement  in  the 
crack  tip  region  as  the  two-dimensional  finite  element  model.  The  lack  of  refinement  in  the  three- 
dimensional  model  requires  that  the  value  of  strain  energy  density  at  the  crack  tip  be  extrapolated 
from  the  strain  energy  density  a  relatively  large  distance  ahead  of  the  crack  tip.  In  order  to  overcome 
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this,  the  strain  energy  density  profile  ahead  of  the  crack  tip  at  the  midplane  of  the  compact  specimen 
is  assumed  to  be  identical  to  die  gradient  determined  from  the  two-dimensional  plane  strain  simula¬ 
tion.  Then,  fracture  initiation  in  the  three-dimensional  simulation  is  determined  to  occur  when  the 
curve  defined  by  three-dimensional  integration  point  data  along  the  specimen  midplane  coincides  with 
the  strain  energy  density  profile  of  the  plane  strain  analysis.  The  computationally  determined  fracture 
initiation  load  is  78.29  kN  (17,600  lb)  and  the  CMOD  is  1.18  mm  (0.0465  in).  The  experimentally 
determined  values  of  fracture  initiation  load  and  CMOD  are  78.06  kN  (17,550  lb)  and  1.17  mm 
(0.0459  in),  respectively. 

PREDICTION  OF  FRACTURE  TOUGHNESS 

Classical  fracture  parameters  are  calculated  from  the  load-displacement  record  and  fracture  ini¬ 
tiation  values  predicted  from  the  finite  element  simulation  using  the  concepts  discussed  previously. 
The  predicted  classical  fracture  parameters  are  compared  with  the  experimental  values. 

The  apparent  stress  intensity  factor,  K,  was  determined  in  accordance  with  ASTM  E-399  specifi¬ 
cations  [18].  While  the  compact  specimen  does  not  allow  for  the  development  of  plane  strain  condi¬ 
tions,  the  specimen  can  be  used  to  determine  the  fracture  toughness  of  22.9  mm  (0.90  in)  thick  sec¬ 
tions  of  HY-100  steel.  The  predicted  K  is  134.1  MPa  (121.7  ksi  'Tin  ).  The  experimental  value  of 
K  is  129.6  MPa  (117.6  ksi  'fin  )  Analytical  and  experimental  values  differ  by  3.5  percent. 

In  a  similar  manner  the  apparent  values  of  J  at  fracture  initiation  were  determined  from  the  two 
load-displacement  curves.  In  accordance  with  ASTM  E-813  procedures  [19],  J  values  are  based  on 
the  load  versus  load  line  displacement  curves.  The  predicted  fracture  initiation  J  value  is  143.7  kJ/m2 
(820  lb-in/in3).  The  experimental  J  value  at  fracture  initiation  is  140.2kJ/m2  (800  lb-in/in3).  Analytical 
and  experimental  values  differ  by  2.5  percent. 

SUMMARY 

Fracture  toughness  at  the  continuum  scale  is  used  as  a  failure  criterion  in  the  current  analysis. 
Material  toughness  at  the  continuum  scale  may  be  appropriately  defined  using  the  strain  energy  per 
unit  mass  at  fracture.  The  continuum  stress-strain  curve  and  the  associated  critical  continuum  strain 
energy  density  are  determined  from  appropriate  analysis  of  observed  deformation  of  simple  material 
tests. 


In  the  current  analyses,  the  continuum  response  used  is  determined  from  cylindrical  tensile 
specimens.  An  accurate  continuum  stress-strain  curve  for  large  strains  which  occur  prior  to  fracture  is 
developed  from  a  hybrid  computational-experimental  procedure.  The  critical  strain  energy  density 
associated  with  fracture  is  determined  from  the  continuum  scale  invariant  stress-strain  curve  and  the 
observed  point  of  failure  from  the  cylindrical  tensile  specimens. 

Nonlinear  finite  element  analysis  of  a  compact  specimen  were  performed  using  the  computer 
code  ABAQUS.  Good  agreement  was  obtained  between  experimental  and  nonlinear  two-dimensional 
finite  element  load-displacement  results  to  approximately  80  percent  of  the  fracture  initiation  load. 
After  this  point,  out  of  plane  displacements  reduce  the  stiffness  of  the  test  specimen.  A  three- 
dimensional  analysis  was  performed.  Nonlinear  three-dimensional  finite  element  results  duplicate  the 
experimental  results  for  the  entire  loading  history. 

Predictions  of  classical  fracture  initiation  toughness  parameters  for  HY-100  steel  agree  to  within 
3.5  percent.  This  demostriates  that  classical  fracture  measures  are  not  fundamental  material  parame¬ 
ters.  Moreover,  these  measures  may  be  derived  through  definition  and  application  of  appropriate  fun¬ 
damental  concepts  and  material  parameters. 
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Two-Dimensional  Finite  Element 


Fig.  6  —  Load-Displacement  Results  Showing  Two-Dimensional  (Plane  Strain)  Finite  Element  Results 


Fig.  9  —  Compact  Specimen  Showing  Curved  Crack  Front 
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Dimensional  Finite  Element  Results 


